Oscillatory dynamics and non-markovian memory in dissipative quantum systems 
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The nonequilibrium dynamics of a small quantum system coupled to a dissipative environment 
is studied. We show that (1) the oscillatory dynamics close to a coherent-to-incoherent transition 
is surprisingly different from the one of the classical damped harmonic oscillator and that (2) non- 
markovian memory plays a prominent role in the time evolution after a quantum quench. 
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To describe dissipation in open quantum systems one 
often relies on phenomcnological approaches. This might 
be sufficient to model experimental observations as long 
as the system is large, the coupling to the environment 
poorly specified, and the accuracy of the measurement 
limited. The rapid progress in controlled access to small 
quantum systems in such distinct fields as condensed 
matter physics, quantum optics, physical chemistry, and 
quantum information science renders a more microscopic 
approach inevitable P, Q • 

Two questions of general current interest are: (1) How 
does the dissipative nonequilibrium dynamics of a quan- 
tum system with only a few degrees of freedom com- 
pare to the standard example of classical dissipation, the 
damped harmonic oscillator (DHO)? Do such systems 
show a coherent-to- incoherent transition of the same type 
as it is found in the classical case? Studying the time- 
evolution of the ohmic spin- boson model (SBM) Q, Q 
we show analytically that the dynamics in the coher- 
ent regime as well as the transition to the incoherent 
one are surprisingly different from their classical counter- 
parts. (2) What is the role of non-markovian terms in the 
dynamics of quantum dissipative systems? By consider- 
ing parameter quenches across the coherent-to- incoherent 
transition in the SBM we show that non-markovian mem- 
ory of the state before the quench heavily affects the time 
evolution after it. When quenching from the incoherent 
to the coherent regime this effect can be so strong that 
the dynamics after the quench is monotonic; the coher- 
ent oscillatory behavior is fully suppressed. We provide a 
qualitative analytical explanation of this numerical find- 
ing. For quenches in the opposite direction nonmonotonic 
behavior can be transferred deep into the incoherent part 
of the dynamics. Recent technological progress to isolate 
small systems from uncontrolled environmental fluctua- 
tions and to manipulate their parameters, e.g. in cold 
atomic gases, led to a considerable interest in quantum 
quenches Q. For the SBM those were not studied so far. 

Model — The SBM is arguably the most important 
model used to describe dissipation in small quantum sys- 
tems coupled to an environment beyond phenomenolog- 



ical approaches. Its Hamiltonian reads 
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with the Pauli matrices cr„, v = x,z and bosonic lad- 
der operators 6? . A spin- 1/2 with Zeeman splitting e 
and tunneling A > between the two states is coupled 
by Afc to a reservoir of bosonic modes with dispersion 
Wfc. The spin-boson coupling is characterized by a spec- 
tral density J(uS) = J2k — gi vcn by the de- 
tails of the microscopic model underlying the SBM [l[. 
The dynamics was mainly studied in the ohmic case with 
J(w) = 2aujQ(uj c - w) and coupling a > Q, 0, i-Q. 
To investigate the fundamental questions raised above 
we reduce the number of parameters by considering the 
(symmetric) case e = 0, the scaling limit, with the high- 
energy cutoff uj c being much larger than any other scale, 
and temperature T = 0. In this strong-coupling limit 
the many-body physics becomes most intriguing and any 
perturbative approach fails [3, Q • It is established that 
the time evolution of the spin expectation value (lp| 
Pit) = (cr z (t)) changes from being coherent, that is 
damped oscillatory, for a < 1/2 to incoherent, that is 
monotonic, for 1/2 < a < 1; at a = 1 a quantum phase 
transition to a localized regime occurs which we do not 
consider. We show that the coherent regime as well as the 
transition still offer surprises. Other parameter regimes 
are investigated in Refs. 

In studies of the nonequilibrium dynamics the system 
is often assumed to be in an initial product state of spin- 
up and the boson vacuum. The time evolution is per- 
formed for a > and P(t) asymptotically tends to zero. 
In the following we refer to this setup as the relaxation 
protocol. We furthermore consider two quench protocols. 
In the first one the system is prepared in the above initial 
state and the time evolution is performed with a coupling 
2«j — 1 > (incoherent regime) up to t — t q . At t q it 
is abruptly switched to otf , with 2a f — 1 = 1 — 2a t < 
(coherent regime) , and the evolution is continued. In the 
second quench protocol we proceed in opposite order and 
quench from 2cti — l<0to2a/ — 1 = 1 — 2ct; > 0. 

Relaxation protocol — The dynamics of the SBM for 
couplings sufficiently away from the transition at a = 1/2 



2 



1(T 



10"' 



10" 



10"' 



10 



10"' 





^a=0.5318 


L a=0.4 






a=0.5~"\ 







10 f a=0.54/ 



no 




fy*a=0.4523 



10 



10 



10" 



• FRG: num. 

— RTRG: num. 

— RTRG: ana. 



10 



FIG. 1: (Color online) Time evolution of a spin-1/2 pre- 
pared in its up state and coupled to an ohmic boson bath 
at t = (spin-boson model). The spin expectation value 
|P(t)| is shown for different couplings a. Dips correspond to 
zeros; to make those visible we use a log y-axis scale. Main 
panel: comparison of the numerical solutions of the FRG and 
RTRG equations as well as the analytical result Eqs. (J2JI- 
Q for a — 0.4682. Right inset: the same as in the main 
panel (various a; analytical expression for a — 0.5318 only for 
which the curves are barely distinguishable). The dynamics 
in the coherent regime is very different to that of the classical 
damped oscillator. Left inset: FRG data for 1 — P(t) at short 
times (dots) compared to the NIB A prediction (dashed lines). 



and outside the transient regime, that is for tA ^ 1. is 
well understood. For small a, that is deep in the coherent 
regime, it compares well to the one of the classical DHO 
[7H9|, [l7[ which is given as the sum of terms each being 
the product of an exponentially damped and an oscilla- 
tory factor. Deep in the incoherent regime the dynamics 
of the SBM is dominated by a single exponentially de- 
caying term 0, [§] in resemblance to that of the DHO. 

We now address if the analogy also holds close to the 
coherent-to-incoherent transition approaching it from the 
coherent side. For a controlled access to this intruiging 
regime we employ the smallness of g = 1 — 2a. Close to 
the transition the characteristic low-energy scale Tk 



A (A/we) ^ 1-0 ' 1 emerges [3, 0). For now we take Tk as 
our unit of energy and measure time in units of T^ 1 . Our 
analytical result for P(t) 
t si 1 reads 



P hc (t)+P p (t) at \g\ < 1 and 
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and the exponential integral E\ . It constitutes one of the 
main results of this Letter. As outlined below it is ob- 



tained employing a two-step procedure which is based on 
complementary renormalization group (RG) approaches: 
the functional (F) and the real-time (RT) [2l| RG. 
In vast contrast to the dynamics of the classical DHO 
which on the coherent side shows infinitely many zeros 
even very close to the coherent-to-incoherent transition 
P(t) only features a single zero for a close to 1/2. This is 
exemplified in the main panel of Fig.[T]which shows |P(t)| 
for a — 0.4682 (solid line with symbols); no additional 
zeros are found for times larger than the ones shown. 
Our analytical result is compared to the numerical solu- 
tions of the differential flow equations (see below) of the 
FRG (symbols) and the RTRG (solid line). The excel- 
lent agreement of all three results indicates that we fully 
control our approximations (see below). The right inset 
shows |P(t)| for different a. In the incoherent regime 
the three curves coincide even better (a = 0.5318). For 
a deeper in the coherent regime further zeros appear, 
e.g. in total two for a — 0.4. 

A 'pole contribution' P p Eq. to P is also found 
in two alternative analytical approaches: the nonintcr- 
acting blip approximation (NIBA) jl[ HI an d conformal 
field theory (CFT) 0. Our ratio n/T{ Eqjffl coincides 
with the one obtained by these methods [22|. While no 
'branch-cut contribution' Pbc Eq. ^ appears in CFT Q 
the one of NIBA leads to a purely algebraic decay which 
turned out to be an artifact of this approximation [2[ . In 
Ref. NIBA was improved leading to 



P hc (t) = -g[l + 3Q(-g)]^. 



(5) 



This result agrees to our Eq. evaluated for t^ 9 ^ ^ 1 
and taking into account that to leading order r| ~ 1 
(that is Tk)- We emphasize that for a < 1/2 the sum 
of the pole contribution Eq. (j3]) and the asymptotic re- 
sult Eq. ([5]) docs not give a meaningful approximation to 
the dynamics for the times shown in Fig. [1] for which the 
coherent and incoherent parts are comparable and non- 
monotonic ('oscillatory') behavior is found. For those it 
is inevitable to keep the branch-cut term in the form of 
Eq. <j2j) . Thus none of the existing analytical approaches 
to the SBM allows to uncover the crucial difference be- 
tween the dissipativc dynamics of the classical DHO and 
the SBM close to the coherent-to-incoherent transition. 
Furthermore, none of the numerical methods applied to 
the SBM [7H9l Il7| was used to investigate the dynamics 
in the transition region. Most likely the data suffer from 
numerical noise which prevents to reliably answer subtle 
questions such as the one of the number of zeros; note 
that the y-axis of Fig. [1] covers six orders of magnitude. 

Further confidence in our RG methods can be gained 
from consideringthe short time dynamics t J$ 1. For this 
NIBA predicts || 1 - P(t) = t 1+9 /T(2 + g) + 0{t 2+2 a) 
which favorably compares to our numerical results (left 
inset of Fig. [TJ. In fact, within both our RG approaches 
this expression can be derived analytically as well. 
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Methods — We next briefly describe our methods. 
Readers interested in results only can skip this part and 
still follow the discussion of the quench dynamics below. 

Using our RG approaches we do not directly study 
the SBM but employ the mapping to the fermionic in- 
teracting resonant level model (IRLM) P, . It is suf- 
ficient to know that the coupling a of the SBM is re- 
lated to the two-particle interaction U of the IRLM: 
1 — 2a = 2U — U 2 = g. The case a = 1/2 corresponds 
to the noninteracting resonant level model (U = 0) and 
is exactly solvable. The FRG flow equations for the 
Kcldysh components of the many-body self-energy of the 
IRLM derived in Refs. HH (relaxation protocol) and [3 
(quench protocols) are controlled to leading order in [/, 
that is g, and can directly be applied to the SBM. 

In the complementary RTRG [2l[ one focuses on the 
reduced density matrix of the local system and describes 
its dynamics in Liouville-Laplace space by 



oo+i0+ 
oo+i0+ 



dEe- lEt U^E), 



(6) 



where H±(E) = [E + iT^E)] -1 is an effective propagator; 
in the IRLM Ti(E) denotes the charge relaxation rate. 
The Laplace variable E can be used as the flow parameter 
25j . leading to 



dT 1/2 (E) 



dE 



= -gT^Uy^E), 



(7) 



with U 2 (E) = [E + iT 2 (E)/2]~ 1 and initial values 
T n {iw c ) = A 2 /u c ; in the IRLM L 2 /2 describes the level 
broadening. Solving Eq. (0 offers the unique possibil- 
ity to identify the singularities of the propagator (poles 
and branch cuts) in the lower half of the complex plane, 
from which the individual terms of the time evolution 
can be identified. The RG equations (0 contain all terms 



0[u£s). However, some terms O 



U 



glectcd. This must be contrasted to our FRG approach, 
which contains all orders in A 2 /oj c . To verify that for 
small \U\, that is small \g\, considering Eq. (0 is suffi- 
cient, we always compare the numerical solutions of the 
FRG and RTRG flow equations (see Figs. HE]). 

In an analytical solution of our RTRG equations (0 
one has to separately consider the branch-cut contribu- 
tion to the integral Eq. ((6]) for all g and the pole contri- 
butions for g > 0. The former follows from the branch 
cut of Ti (E) on the imaginary axis starting at the pole 
-iT* 2 /2 of U 2 (E). Taking T 2 {E) « T 2 in IL 2 (E), we get 
ri(£) w (T* 2 /2 - iE)-9 from Eq. (0, leading to 



Pbc(t) 



e 2 



V' 2 t/2 f oo 



dx Im 



-f +x- e n s x -9 
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Replacing x g by t 9 in the denominator gives a very 
good approximation for t ^ 1 and leads to Eq. (0. A 
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FIG. 2: (Color online) Spin expectation value P(t) of a 
spin- 1/2 with the initial time evolution given by the spin- 
boson Hamiltonian and a coupling from the incoherent regime 
oti > 1/2. At t = t q it is switched to a/ < 1/2 from the co- 
herent one. The scaling limit is realized by A/oj c — 1/200. 
Times are restricted to t > t q ; for < t < t q , see Fig. [1] Left 
inset: comparison of numerical FRG and RTRG data; the 
non-markovian memory completely suppresses the coherent 
'oscillatory' behavior. It only survives for ctj smaller than a 
critical coupling a c ; see the upper [numerical solution of FRG 
equations], central [numerical solution of RTRG equations], 
and lower [analytical result Eq. (|10p ] panels. Right inset: de- 
pendence of a c on the time T l K t q the evolution was performed 
in the incoherent regime (FRG data). Two different defini- 
tions of a c are compared. In the first a c is defined as the q/ 
at which the first zero of P(t) can be observed, in the second 
as the a/ at which the curve first shows a nonmonotonicity. 
The two values barely differ. 



comparison with the numerical solution shows that to 
find the pole positions z\ = ±fi — iT'l (z 2 = —iT 2 /2) 
of ILi(E) (a 2 (E)) one can set T^E) w T 2 (E) iz 1 
(TxiE) » T 2 (E) « T*) in U n (E). The RG equations can 
then be solved analytically leading to Eq. (0 as well as 
the pole contribution to P(t) Eq. (0. 

Quench protocol 1 — We next investigate the role of 
non-markovian memory in dissipative dynamics. In 
Fig. [5] we show P(t) when quenching at t q from the in- 
coherent to the coherent regime. The two different cou- 
plings on and af lead to the two characteristic scales 
T % K and Tj£, which in the scaling limit differ by orders 
of magnitude. Therefore the model parameters A and 
lo c cannot be scaled out as efficiently as in the relax- 
ation protocol (by taking Tk as the unit of energy). To 
minimize the influence of the transient dynamics before 
the quench we consider T l K t q between 5 and 10. In the 
left inset we compare the data obtained by the numer- 
ical solution of the FRG [24[ and RTRG equations (sec 
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below). The agreement for cvj/j close to 1/2 is excel- 
lent. Even though the same coupling a/ = 0.4682 as in 
the main panel of Fig. Q] is considered the data for the 
quench dynamics displayed in the left inset of Fig. [5] are 
monotonic; no indication of coherent behavior is found. 
The incoherent dynamics before the quench thus heavily 
affects the one afterwards. To further investigate this, 
results for smaller af are shown in the upper (FRG) and 
central (RTRG) panel. Nonmonotonic ('oscillatory') be- 
havior can only be found for a / being smaller than some 
critical coupling a c . As |1 — 2a c \ <C 1 docs not hold 
strictly, our two approximate RG methods give values 
for a c which differ by a few percent. Its precise value 
can only be obtained using a method which treats higher 
order contributions consistently 

Quenches were so far not studied using RTRG; the 
technical details will be given elsewhere. In this approach 
the memory of the spin state is preserved in the exci- 
tations of the bath. Classifying the memory contribu- 
tions by the number of bath excitations that have gone 
through the quench one can show that the terms with 
n excitations are proportional to A" = {T l K lT^) n / 2 
(A/w c ) 2ra ' 9/ !, with A<1. Restricting ourselves to the 
case of a single excitation with Matsubara frequency A 
we obtain {t 1 = t — t q > 0) 

/>oo 

P(t')=pf(t')P i (t q )-9i dAF[(t')Fi(t q ), (9) 
Jo 

with k = i, f and P K , T^ 2 computed as in the relax- 
ation protocol with the corresponding on if and initial 
condition P K (0) = 1. The first term describes standard 
relaxation while the second memory term implies (dissi- 
pative) non-markovian dynamics. It has no analog for 
quenches in closed systems. The RTRG data of Fig. [2] 
were obtained by numerically performing the E and A 
integrals in Eq. (|9|) on the basis of the numerical solu- 
tion of Eq. ([7]). Fixing uj c the results depend only very 
weakly on A 2 /w c (via A). Varying (A/w c ) 2 by an order 
of magnitude around the average value 2.5 • 10 -5 results 
in a change of a c by only a few percent. 

To gain qualitative analytical insight of the interplay 
of the relaxation dynamics and the non-markovian cor- 
rection we evaluate the integrals keeping only the terms 
which dominate for \gi//\ <C 1. The result 

W- « 2(1 - A)e- r ' tt ' cos(^t') + Ae- r i ft '/ 2 (10) 

is shown in the lower panel of Fig. [2] The memory gener- 
ates a coherent and an incoherent contribution oc A. The 
first has a negative sign and suppresses the coherent part 
while the second enhances the incoherent term compared 
to the one of the relaxation protocol; for <?/ <C 1 the lat- 
ter is subdominant and thus not written in Eq. (|10p . This 
explains the appearance of a critical a c . 




FIG. 3: (Color online) The same as in Fig. [5] but for quenches 
from the coherent to the incoherent regime. The spin expec- 
tation value |P(4)| for on = 0.4682 and different quench times 
t q (indicated by the vertical dotted lines) is shown. At the 
respective t q the i-axis scale is switched from T^t to T K t. 
Main panel: data of the numerical solutions of the FRG and 
RTRG flow equations. The thin dashed line is an exponential 
term with rate and a subleading power-law correction. 

Inset: FRG data with and without the memory term. The 
non-markovian memory leads to a transfer of nonmonotonic 
behavior deep into the incoherent part of the dynamics. 



Quench protocol 2 — We finally discuss the dynamics 
when quenching at t q in the opposite direction, that is 
from the coherent to the incoherent regime. To keep the 
discussion transparent we focus on o^'s for which only a 
single zero at time to is found in the relaxation protocol. 
In Fig. [3] we show \P(t)\ obtained from the numerical 
solution of FRG and RTRG equations for a = 0.4682. 
In the relaxation protocol TkU) w 5 for this coupling 
(see Fig.[T]). Different quenching times t q are considered. 
For t q > to the dynamics at times larger than t q very 
quickly adapts to the new rate « of the incoherent 

dynamics [see Eq. ([5|) and the dashed line in Fig. |3] . The 
behavior is significantly different for t q < to- In this case 
nonmonotonic ('oscillatory') behavior is found for times 
t > t q at which the time evolution is performed with 
a f = 0.5318 > 1/2. The smaller t q - t the further the 
memory term (see the inset of Fig. [3|) transfers the zero 
into the incoherent regime. 

Summary — We have shown that the coherent dynam- 
ics of the ohmic SBM, being the prototype model of 
dissipativc quantum mechanics, close to the coherent- 
to-incohcrcnt transition is very different from the one 
of the classical DHO. Our study furthermore revealed 
the crucial importance of non-markovian memory in the 
noncquilibrium time evolution when quenching across 
this transition. Both these findings come as surprises. 
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